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Abstract 

Let us denote by p(x|K) the space density of the points where 
identical particles of some kind, e.g. vr^ mesons, with momentum K 
are produced. When using the HBT method to determine p(x|K) 
one encounters ambiguities. We show that these ambiguities do not 
affect the even cumulants of the distribution p(x|K). In particular, the 
HBT radii of the homogeneity regions, which are given by the second 
order cumulants, and the distribution of distances between the pairs 
of production points for particles with momentum K can be reliably 
measured. The odd cumulants are ambiguous. They are, however, 
correlated. In particular, when the average position (x)(K) is known 
as a function of K there is no further ambiguity. 
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1 Introduction 



Femtoscopy, i.e. the study of momentum correlations among particles at 
small relative velocities in order to get information about the interaction re- 
gions where the hadrons are produced, is now considered so important that 
workshops devoted to correlations and femtoscopy are organized every year 
(see e.g. [1]). The HBT method is the oldest and much used tool of fem- 
toscopy. For reviews with hundreds of references see e.g. [2l[3llll[5]. There is 
a variety of models which predict the correlation functions depending on some 
parameters. Fitting these parameters to the data one obtains information 
about the interaction regions. Much less attention has been devoted to the 
problem: given the data what can be deduced without using a specific model. 
A well-known, important result is that very little can be learned about the 
whole interaction region. It is much more fruitful to concentrate on the ho- 
mogeneity regions. Homogeneity region K is the region where hadrons with 
momentum K are produced. For a homogeneity region, much information 
with little model dependence can be obtained using the imaging method (see 
[6j and references quoted there). The imaging method, however, yields the 
distribution of distances between the points where identical hadrons with a 
given momentum are produced. This is less information than contained in 
the distribution of the production points themselves. Therefore, the question 
arises: can one go any further and if not what information from models must 
be added? 

Let us denote p(x|K) the probability distribution for the points where 
the hadrons (of the type considered) with momentum K are produced in 
a given kind of interaction. We will call p(x|K) the profile of homogeneity 
region K. Barring some pathological cases, the profile can be unambiguously 
determined when all its moments are known. All the moments can be un- 
ambiguously determined when all the cumulants are known. In the present 
paper we concentrate on cumulants. We find that all the even cumulants 
can be unambiguously determined from the data - i.e they are measurable. 
This implies in particular two important, known results. The famous HBT 
radii and the distribution of distances between the production points are 
measurable. Actually, knowing the distribution of distances is equivalent 
to knowing all the even culmulants. A single odd cumulant is not measur- 
able. In particular, the K dependence of the first order cumulants, which are 
equal to the components of the the first order moments (x)(K), can have an 
arbitrary dependence on K. This implies, as a very special case, the well- 
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known result that the position of the center of the whole interaction region is 
not measurable. Changing the relative positions of the homogeneity regions, 
however, one usually has to change their shapes accordingly. Unless p(x|K) 
is a Gaussian, rigid shifts of the homogeneity regions with respect to each 
other are not allowed. For an example where the skewness gets changed see 
[7|. The common knowledge is (cf. e.g. [^) that the correlation functions 
can provide at best the distribution of the relative positions of particles with 
identical velocities. We find that for given (x)(K) the full profiles p(x|K) 
can be measured. This result may be useful. Macroscopic models, like the 
hydrodynamical ones, are more likely to predict correctly the K-dependence 
of the positions of the homogeneity regions than finer details of their shapes. 
According to our result, however, function p(x|K) is measurable when the 
input includes (x)(K) for every K. 

Our analysis is based on two main assumptions. One is that formula 
is approximately valid. This can be interpreted as the assumption that, 
besides the Bose-Einstein correlations, the particles are uncorrelated. The 
other assumption is that formulae ([2]) and Q are approximately consistent. 
This is known to be true when we are sufficiently close to the classical limit. 
Both these assumptions hold for a broad class of models. Our results are 
model independent for models within this class. 

2 Basic formulae and assumptions 

The probability of emitting at time t a particle with momentum K from 
space-time point X is given by the emission function S{X,K) [3], fS]. Since 
the particles are emitted on mass shell, we have 



where m is the particle mass. With this definition, the unnormalized profile 
of homogeneity region K, i.e. the space distribution of points where the 
identical hadrons, say tt~^ mesons, with momentum K have been produced, 
is given by the integral over time 



This profile differs only by an X- independent factor from the profile p(X|K), 
normalized to one at each K. Therefore, either profile can be used to describe 
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the homogeneity regions. In this section we changed the notation for the po- 
sition vector from x to X. This has been done to improve the correspondence 
with the other definition of the emission function described below. Actually, 
unless the sources are point-like in space-time, X in the following is not quite 
the position vector of a particle, but our basic assumption is that the two 
definitions of S are, with sufficient accuracy, consistent. Then, except in the 
integrands of next three formulae, x and X can be used exchangeably. 
On the other hand one uses [1], [8] the formula 

S{X,K) = J2 [d'yT*{X + ^y)TF{X-^-y)e'^y, (3) 

F 

where the sum is over the states of all the other particles in the system. 

Following [8] we make the approximation that each source can emit a par- 
ticle only at one time, though different sources may emit particles at different 
times. This assumption of instant sources cannot be rigorously true, it would 
contradict Heisenberg's uncertainty relations, but as an approximation it is 
in the spirit of the models being used. Thus e.g. Kopylov and Podgoretsky in 
their classical paper [9] explicitly introduce long-hved sources, but then they 
average over an internal, unmeasurable parameter [E) which makes these 
sources equivalent (for the HBT analysis!) to sets of instant sources. Many 
papers go even further and use sources which are, or can be decomposed 
into sources, point-like in space-time [10\. A general derivation of the 
relation between correlation functions and emission functions, which is still 
[B] considered to be the most detailed, has been given by Anchishkin, Heinz 
and Renk [TT]. Analyzing it one easily finds that the proof breaks down if 
the sources are not, or at least for the HBT analysis cannot be replaced by, 
instant sources. One also finds [12j that S for coherent sources extended in 
time cannot be interpreted as a phase space density of the particles produced 
at a given time. 

For instant sources the time component of y is zero by assumption. 
Therefore, the integration over y'^ and the term K^y'^ in the exponent drop 
out. Denoting by ti the time when a source fires and using the labels Fi for 
the sources that fire at t = U we get 

S{X,U,K) = [d'y ^T;^(X + iy)T^,(X-ly)e-^-y. (4) 

F 

Usually ti is considered a continuous variable, but for comparison with the 
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Wigner function it is more convenient to start with discrete time variables 
and, if necessary, introduce the continuous variables later. 

The sum over Fi yields a result proportional to the complex conjugate of 
the corresponding density matrix. Denoting the (real, non-negative) propor- 
tionality coefficients by N{ti) one gets 



S{y.,U,K) = Nit,) jd'y p%X + ly,X~ ^y))e-^K-y = N{ti)W{X,U,K). 

(5) 

Actually, the integration d^y yields the complex conjugate of the Wigner 
function. Since, however, the Wigner function is real the second equality is 
correct. Summation over tj yields 

p(X|K) = iViy(X,K), (6) 

where the Wigner function W(K, K) is the average over the Wigner functions 
W(K, ti, K) and is a K-dependent normalization factor chosen so that for 
every K the probability distribution p(X|K) is normalized to one as it should. 

This result shows that the definition of the emission function as a phase 
space density, used to derive formula (I2l), is only approximately consistent 
with ([3]). A Wigner function is not a classical distribution. It can take nega- 
tive values, it is bounded by the condition K)| < n^^ etc. Moreover, 
condition ([1]) is automatically satisfied for the density, while for the Wigner 
function = ^(p? + P2), where both pi and p2 are on mass shell, can sat- 
isfy ([T]) only approximately for small |q|. Using the Wigner function for the 
phase space distribution can be justified in the classical limit [13] • Crude 
estimates (cf. e.g. [H]) suggest that it is a very good approximation for the 
homogeneity regions in heavy ion scattering, though not necessarily so for 
much smaller interaction regions as e.g. in pp or e~^e~ scattering. 

The standard relation between Wigner functions and density matrices 
yields 

p(K,q) = p(K,0) Jd'Xpi^\K)e~'^ '', (7) 

where 

K = ^(Pi + P2), q = Pi-P2 (8) 
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and p(K, q) is the density matrix corresponding to the Wigner function 
Vr(X,K). 

Using the HBT method one gets information about the density matrix 
p(K, q) from the measured two-, and sometimes more-, particle correlation 
functions. In most models, when final state interactions including resonance 
production are neglected or corrected for, the basic formula is [T3] 

n 

P j=l 

where P(pi,...,p„) is the n-particle probability distribution, the summa- 
tion is over the n\ permutations of the indices j and C„ are normalization 
constants. 

When this density matrix is known, the profiles of the homogeneity re- 
gions p(X|K) can be obtained by inverting the Fourier transformation. There 
are, however, ambiguities in the determination of the density matrix from the 
data. The fits to the data do not change under the transformations [7], [16] 

p(K, q) ^ e^^(PiV(K, q)e-'^^P^\ (10) 

where / can be any real-valued function. They are also invariant under 
complex conjugation of the density matrix, but this is less interesting since it 
corresponds just to a space inversion of the homogeneity region. Conversely, 
(fTOj) and complex conjugation are the only transformations of the density 
matrix which leave the probabilities (jH]) unchanged [TGJ. The purpose of the 
present paper is to quantify the ambiguities in the profiles p(K\K.) due to 
the ambiguity (fTU]l . 



3 Results 

In the terminology of the calculus of probabilities, relation ([7]) means that, 
at given K, p(K, q) / p(K, 0) is the characteristic function of the distribu- 
tion p(x|K). The expansion of the logarithm of the characteristic function 
in powers of the components of q yields the cumulants }C{rx,ry,rz) of the 
distribution p(x|K): 

logp(K,q)= ^^%/C(r.,r„r,)+logp(K,0). (11) 

^.T* ^11' ^Z* 
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We introduce the notation r = + Ty + and call a cumulant even (odd) 
when r is even (odd). Let us rewrite the logarithm of the density matrix in 
the form 

logp(K, q) = logpo(K, q) + «x(K, q), (12) 

where both x and po are real. The hermiticity of the density matrix implies 
that po is an even function of q and x is an odd function of q. Thus, po 
yields the even cumulants and x the odd cumulants. Since the ambiguity 
f fTOj) affects only x, the even cumulants can be unambiguously found using the 
HBT method. This is our first result. A notable case are the r = 2 cumulants 
which yield the famous HBT radii (see e.g. [2]). The measurability of the 
HBT radii has been already pointed out in [7j. Another useful implication 
is the measurability the distribution of distances between pairs of points 
where particles with given momentum K are produced. The corresponding 
characteristic function is 



0(q|K) = |d3Xid3x,p(Xi|K)p(X,|K)e-(^-^^) = ^^^^^^^^y""^ ■ 

The log0(q|K) is an even function of q, therefore, all the odd cumulants van- 
ish, while the even ones are just twice the measurable cumulants of p(x|K). 
This well-known result is basic for the method of imaging (cf. e.g. ^6j). In 
the present derivation the smoothness assumption is replaced by the closely 
related assumption that formulae ([2]) and ^ can be simultaneously used. 

The general formula for the unmeasurable contributions from function / 
to the cumulants is 



SlCir r r)=(—\ ^ ^^^^ ^ " ^'^^ (U) 



'-tV-' y/(K) !-(-! 
2 J OKl-dKl'dKl^ 

The r = 1 cumulants give (x)(K), i.e. the position of the center of the 
homogeneity region as a function of K. The unmeasurable parts of these 
cumulants yields for the uncertainty of position 

5(x)(K) = V/(K). (15) 

Formally, this is analogous to the uncertainty of the electromagnetic vec- 
tor potential due to the freedom of time-independent gauge transformations. 



7 



E.g. it is possible to choose /(K) so that for every homogeneity region 
{z) = 0. The difference is, however, that in electrodynamics the choice of 
gauge does not affect the physics and is purely a matter of convenience, while 
here only one choice of /(K), except for an irrelevant additive constant, is 
correct and yields the true profiles of the homogeneity regions. Since ex- 
perimental momentum distributions give no information about /(K), this 
function must be deduced from theory or from a model. As seen from f[T^ . 
it is enough to fix the function (x)(K), because then 5(x)(K) = 0, / reduces 
to a constant and all the cumulants can be unambiguously calculated. Con- 
sequently, the full profiles p(x|K) are measurable. It should be kept in mind 
that changing /(K) does not mean only shifting the homogeneity regions 
with respect to each other. In general, the higher cumulants also change 
and consequently the homogeneity regions are not only shifted, but also de- 
formed. A one-dimensional example where the change affects the first and 
third order cumulants has been discussed and illustrated by a figure in [7]. 

From the two-body correlation function one can obtain po(K,q) - the 
absolute value of the density matrix. Often it is assumed, usually implicitly, 
that po can be used as an approximation for the true density matrix p. As 
seen from our analysis, po yields correctly all the even cumulants and thus 
embodies all the information contained in the distribution of distances be- 
tween the points where the hadrons have been produced. Replacing p by po 
implies also, however, putting arbitrarily all the odd cumulants equal zero. 
This may lead to disagreement with experiment for the three-particle cor- 
relation function (cf. e.g. the review [2J ). If the three-particle correlation 
function happens to be reproduced correctly and if assumption Qj holds, 
there will be also agreement with experiment for all the more-particle cor- 
relation functions The profiles may be still wrong, however, because of 
the ambiguity ffTU]) . 

4 Conclusions 

Let us summarize our conclusions. All the even cumulants of the profiles 
p(x|K) are unambiguously measurable. This implies in particular the mea- 
surability of the HBT radii and of the distributions of distances between 
points where the identical particles with the same momenta are produced. 
These results are known, but we give a new, very simple proof which clearly 
exhibits the necessary assumptions. 
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The odd cumulants, though separately not measurable, are strongly cor- 
related with each other. In particular, it is enough to know the positions of 
the centers of all the homogeneity regions (x)(K), in order to reconstruct the 
full profiles p(x|K). Our analysis justifies the following recipe. 

Prepare as input: the data on the single-particle momentum distribu- 
tions, the two-body and the three-body correlation functions and a distribu- 
tion (x)(K) assumed to be correct. More-particle correlation functions could 
be used to reduce the statistical errors, but in principle they are redundant. 
Use any model which reproduces correctly the experimental input to find 
a density matrix pm(K,q). This will be one of the infinitely many density 
matrices giving exactly the same fit to the data, but in general very different 
profiles of the homogeneity regions. Use this matrix to calculate the three 
components of the first moment (x)j„(K). The subscript in is used to dis- 
tinguish the moments given by the model from the true moments contained 
in the input. Calculate the correction 

5(x)(K) = (x)(K)-(x),„(K) (16) 

Substitute it into (ITBi) and solve for /(K). Correct the density matrix pj„ 
according to ffTOj) . Substitute the corrected density matrix into ([7]) and solve 
for p(x|K) by inverting the Fourier transformation. 
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